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Recent numerical simulations with different techniques have all suggested the existence of a contin- 
uous quantum phase transition between the Z2 topological spin liquid phase and a conventional Neel 
order. Motivated by these numerical progresses, we propose a candidate theory for such ^2— Neel 
transition. We first argue on general grounds that, for a SU(2) invariant system, this transition 
cannot be interpreted as the condensation of spinous in the Z2 spin liquid phase. Then we propose 
that such Z2— Neel transition is driven by proliferating the bound state of the bosonic spinon and 
vison excitation of the Z2 spin liquid, i.e. the so called (e,m)— type excitation. Universal critical 
exponents associated with this exotic transition are computed using 1/7V expansion. This theory 
predicts that at the Z2— Neel transition, there is an emergent quasi long range power law correlation 
of columnar valence bond solid order parameter. 
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I. INTRODUCTION 

Thanks to the rapid development of numerical tech- 
niques, more and more candidates of exotic liquid states 
have been identified in frustrated spin models^"— , hard- 
core quantum boson model^i^, or even Hubbard modelii. 
All these phases that are identified numerically are fully 
gapped liquid phases with short range correlation be- 
tween both spin order parameters and also valence bond 
solid (VBS) order parameters. The simplest fully gapped 
spin liquid state is the Z2 topological liquid state, which 
has the same topological order as the toric code model—. 
In addition to the fully gapped spectrum, the computa- 
tion of critical exponents at the order-disorder transition 
of these model o'^'^'^" , and the computation of topological 
entanglement entropy^ both convincingly proved that 
the spin liquid states of some of these models (such as 
the Ji — J2 model on the square lattice, and the extended 
Bose-Hubbard model on the Kagome lattice) are indeed 
the Z2 topological liquid. In some other models (such 
as the spin- 1/2 Heisenberg model on the Kagome lat- 
ticeii^, and the Hubbard model on the honeycomb lat- 
ticeii), although an accurate topological entanglement 
entropy computation is still demanded, it is broadly be- 
lieved that the spin liquid state is indeed the Z2 liquid 
state, or a similar (Z2)" liquid state. 

Besides the spin liquid state itself, the quantum phase 
transitions of these models are equally interesting. For 
example, continuous quantum phase transitions between 
Neel order and a fully gapped spin liquid phase have been 
found in the honeycomb lattice Hubbard modelii, and 
the Ji — J2 spin-1/2 Heisenberg model on the square lat- 
tice^i^. In terms of the Landau- Ginzburg (LG) theory, 
this transition should be an ordinary 0(3) transition, and 
the Z2 liquid phase is identified as the disordered phase, 
while the Neel phase is the ordered phase. However, be- 
cause the Z2 liquid phase has a nontrivial topological 
order and topological degeneracj*i^, it cannot be adia- 
batically connected to the trivial direct product state. 



thus it should not be identified as the trivial disordered 
phase in the classical case. Thus if the Z2— Neel transi- 
tion exists, it means that the quantum disordering of the 
Neel order and the emergence of the Z2 topological or- 
der happen simultaneously at one point, this unusual fact 
implies that this Z2— Neel quantum critical point (QCP) 
must be an unconventional one that is beyond the LG 
paradigm. The goal of this paper is to understand this 
unconventional QCP. 



II. FAILURE OF THE SPINON THEORIES 

We first argue on general grounds that such a contin- 
uous ^2— Neel transition cannot be understood using an 
ordinary spinon theory. We stress that we will only con- 
sider SU(2) invariant systems here. 

First of all, if this Z2 spin liquid phase has a 
gapped fermionic spinon excitation /„, then a Neel or- 
der parameter can in principle be represented as A^ ^ 
{~^y fi a'^apfi,p- Thus it appears that we can interpret 
this Z2— Neel transition as the disorder-order transition 
of the vector N using an ordinary Landau-Ginzburg the- 
ory. However, this theory is incorrect because the vector 
N does not carry any gauge charge, thus the order pa- 
rameter does not immediately suppress the Z2 topologi- 
cal order. This implies that between the Z2 spin liquid 
and the Neel order with nonzero (A), there must be an 
intermediate state with the coexistence of both Neel or- 
der and Z2 topological order, and it is usually called the 
Neel* state. Thus a direct continuous transition between 
Z2 and Neel order cannot be obtained this way without 
fine-tuning. 

In order to suppress the Z2 topological order, the usual 
wisdom is to condense a topological excitation that car- 
ries the Z2 gauge charge. Then after the topological ex- 
citations are condensed, the Z2 gauge field is Higgsed, 
and the topological order disappears. Along with sup- 
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pressing the topological order, if we want to induce spin 
order simultaneously, then the excitation that condenses 
must also carry certain representation of the spin SU(2) 
symmetry group, in addition to the Z2 gauge charge. 
Let us call this gauge-charged spin excitation as spinon 
in general. Then the nature of the spin order and the 
universality class of this transition both depend on the 
particular spin representation of spinon. 

The smallest representation of SU(2) is spin-1/2 rep- 
resentation, and there is no consistent "fractional" rep- 
resentation of SU(2) group that is smaller than spin-1/2. 
Thus let us first assume the spinon is a spin-1/2 boson, 
which is described by a two component complex boson 
field Za — (zi,Z2)*, and Za is subject to the constraint 
•Zip + |-Z2p = 1- Then Za is coupled to a Z2 gauge field 
in the following way: 

H^J2J2 + H.c. + ■■■ (1) 

where the ellipsis stands for higher order interaction 
terms. ct|_^ is the Z2 gauge field that is defined on the 
link (i, fi) of the lattice, and Eq. [T]is invariant under the 
gauge transformation 

Zi,a VtZi,a, cr-.^, ■q,(jl^-qi+^,, (2) 

where rji = ±1 is an arbitrary Ising function defined on 
the sites of the lattice. The condensed phase of z^ is 
the spin ordered phase, and because Za is coupled to the 
Z2 gauge field, the Z2 topological order is automatically 
destroyed due to the Higgs mechanism in the conden- 
sate of Za ■ The gapped phase of is the deconfined Z2 
topological phase. 

Since Zq. has in total two complex bosonic fields, i.e. 
four real fields, then with the constraint |zip -I- |z2p = 1, 
the entire configuration of Zq. is equivalent to a three di- 
mensional sphere S^. Since the spinon field Za is coupled 
to a Z2 gauge field, then the physical configuration of 
the condensate of Zq, is S'^ /Z2, which is mathematically 
equivalent to the group manifold S0(3). Since Za itself 
is not a physical observable, inside the condensate of z^ 
the physical observables are the three following vectors: 

Ni = Re[z*icr2'az], N2 = lTa[zH(jy az], N3 = zVz. (3) 

A simple application of the Fierz identity '^ap'^ip ~ 
'^SapS^'y — SajsS-yp proves that these three vectors are or- 
thogonal with each other. Since the first homotopy group 
of SO (3) is TTi [SO (3)] = Z2, inside this spin ordered phase 
there are vortex-like topological defects. Two of these 
vortices can annihilate each other. 

The spin-1/2 boson field z^ can be viewed as the low 
energy mode of the usual Schwingber boson ba, but our 
argument is more general, and it is independent of the 
microscopic origin of Zq . If we identify one of the three 
vectors Ni as the Neel vector, then this phase must have 
two other spin vector orders that are perpendicular to 
the Neel vector. The condensation transition of Zq while 



coupled to a Z2 gauge field is usually called the 0(4)* 
transition^"^. 

Now let us assume the spinon of the Z2 topological 
phase carries a spin-1 representation. A spin-1 represen- 
tation is a vector representation of SU(2), i.e. it can be 
parametrized as a unit real vector n, |np = 1. Now the 
coupling between the spinon and Z2 gauge theory reads 

^ = EE-Km«X+m + --- (4) 

a 

Again, since n couples to a Z2 gauge field, it is not a 
physical observable: n and — n are physically equivalent. 
If vector n condenses, the condensate is in fact a spin 
nematic, or quadrupole order, with physical order pa- 
rameter 

Q-b = n-n' - Uab. (5) 

This spin order has manifold S'^ /Z2, which also supports 
vortex excitation since tti [6*^/^2] — Z2. One example 
state of this type is the spin quadrupolar state that has 
been observed in the spin-1 material NiGa2S4i^. 

We have discussed two types of unconventional QCPs 
between Z2 liquid phase and spin orders. In cither case, 
the spin ordered phase is different from the ordinary 
coUinear Neel order, because a Neel order should have 
ground state manifold (GSM) S^. In particular, in both 
cases we have considered, the spin ordered phase must 
have a nontrivial homotopy group tti , which corresponds 
to the vison excitation of the Z2 gauge field. Gener- 
alization of our analysis to higher spin representations 
is straightforward, but the conclusion is unchanged. As 
we already discussed, in Ref.-^ and Ref.^"'^, a continuous 
quantum phase transition between a fully gapped spin 
liquid phase and a Neel order was reported. If the fully 
gapped spin liquid discovered in these numerical works 
is indeed a Z2 spin liquid as we expected, then such con- 
tinuous quantum phase transition is beyond the spinon 
theory discussed in this section. In order to understand 
the continuous transition between the gapped spin liq- 
uid and Neel order reported in the phase diagram of the 
Hubbard model on the honeycomb latticeii, in Ref J^"— 
the authors had to introduce extra "hidden" order pa- 
rameters in the Neel phase, which change the GSM of 
the Neel phase completely. 

In this section we argued on general grounds that the 
^2— Neel transition cannot be interpreted as the conden- 
sation of an ordinary spinon. Our argument is indepen- 
dent of specific spin model or lattice structure. However, 
this argument can only be applied to SU(2) invariant 
systems. For a system with U(l) symmetry, for exam- 
ple the hard-core Boson model on the Kagome lattice 
discussed in Refi^'Si'' 1^, the transition between Z2 topo- 
logical phase and the superfluid phase can be understood 
as the condensation of a fractionalized "half-boson" that 
couples to the Z2 gauge field, and this transition is the 
so-called 3d XY* transition. 
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III. EXOTIC Za-NEEL QUANTUM CRITICAL 
POINT 

A. Phase diagram around Z2 spin liquid driven by 
e and m excitations 

In order to understand the direct continuous transi- 
tion between the Zi spin liquid and the Neel phase, we 
should first put these two phases in the same phase di- 
agram. One candidate theory that contains both phases 
was proposed in Ref.^*^. Let us first write down a minimal 
unified field theory proposed in Ref.— : 

Q — 1 a — 1 

+ Sz\Za'^ ^ Sv\Va'^ + -^u.vpau^vbp^ (6) 

TT 

In this field theory, there are two types of matter fields, Za 
and Va , and they are interacting with each other through 
a mutual Chern-Simons (CS) theory, which grants them 
a mutual semionic statistics i.e. when adiabatically 
encircles through a closed loop, the system wave- 
function acquires a minus sign. This is one of the key 
properties of the topological phase. Here z^ corre- 
sponds to the electric (e— type) excitation of the Z2 liq- 
uid, and Va corresponds to the magnetic (m— type) exci- 
tation. Va is usually called the vison excitation. 

The minimal field theory Eq. |6] has symmetry 
SU(A'2)xSU(7V„). However, depending on the details of 
the microscopic model, the higher order interactions be- 
tween matter fields can break this symmetry down to its 
subgroups. We will first ignore this higher order symme- 
try breaking effects, and focus on the case with iV^ = 2, 
and iVt, = 1. In Ref.™, the authors used the model Eq. [6] 
with = 2, iVi, = 1 to describe the global phase dia- 
gram of spin- 1/2 quantum magnets on a distorted tri- 
angular lattice, which is a very common structure in 
many materials. The same theory can be applied to the 
square and honeycomb lattice as well, and in this paper 
we will take the square lattice as an example. Here z^ is 
a bosonic spin-1/2 spinon, and v is the low energy mode 
of a vison, and it corresponds to the expansion of the 
vison at two opposite momenta ztQ: 

T ^ ve'^-^ -Vv* e-'^-\ (7) 

thus f is a complex scalar field. On the square lattice or 
distorted triangular lattice, there is a Zs anisotropy on 
that is allowed by the symmetry of the lattice^"'^^. This 
anisotropy is highly irrelevant in the quantum critical 
region, and it will be ignore throughout the paper. 

The phase diagram of this model is tuned by two pa- 
rameters: Sz and s„, and depending on the sign of these 
two parameters, there are in total four different phases 
(Fig.©: 

Phase 1. This is the phase with Sz > 0, Si, > 0. In 
this phase, both matter fields Za and v are gapped, and 
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FIG. 1: The global phase diagram of Eq. (6] which describes 
four different states on a distorted triangular lattice, or a 
square lattice. Eq. [6] assumes that the e— type excitation Za 
and m— type excitation v condense separately. 

they have a topological statistic interaction through the 
mutual CS theory. Since all the matter fields are gapped, 
the low energy properties of phase 1 is described by the 
mutual CS theory only. The mutual CS theory defined 
on a torus has a four-fold degenerate ground state, thus 
this phase is precisely the gapped Z2 topological phase^^. 

Phase 2. Sv > 0, Sz < 0. When vison v is gapped, 
integrating out v induces a Maxwell term for gauge field 
6^, which implies that the flux of 6^ is condensed. In 
other words the flux-creation operator (denoted as Aib) 
acquires a nonzero expectation value. Mj, corresponds to 
a Dirac monopole configuration of in the space-time. 
Due to the mutual CS coupling between gauge fields 
and 6p, the condensate of A^b breaks to a Z2 gauge 
field. Thus after we integrate out v and b^, the spinon Za 
is only coupled to a Z2 gauge fields. Thus when Nz = 2, 
the condensate of Za has GSM SO (3) as was discussed 
in the previous section. An example of this phase is the 
spiral spin density wave phase. Once we assume > 0, 
Eq. [S] precisely reduces to the previously studied 0(4)* 
theory for the transition between Z2 spin liquid and spiral 
spin order^'^. 

Phase 3. < 0, Sz > 0. This is a phase where 
V condenses while z^ is gapped out. This phase is the 
four fold degenerate columnar VBS phase that breaks the 
reflection and translation symmetry of the lattice. The 
columnar VBS order parameter can be written as v^Ma, 
where ^Aa is the monopole operator of gauge field a^, 
which creates a 27r flux of a^. When Sz > 0, spinon z^ 
is gapped, and it leads to a Maxwell term for a^. This 
implies that M-a is condensed, and it breaks b^ to a Z2 
gauge field. In this case the low energy effective theory 
that describes phase 3 is a complex field v that couples 
to a Z2 gauge field, thus our theory reduces to the pure 



4 



vison theory that was thoroughly discussed in Ref<2i. 

Phase 4- < 0, Sz < 0. This is a phase where 
both Za and v condense. Because in this phase the only 
gauge invariant order parameter that condenses is A'^ ~ 
z'^ffz, this is precisely the collinear Neel phase with GSM 
S^. In fact, when v is condensed, the gauge field 
acquires a mass term 6^ due to the Higgs mechanism. 
Then integrating out v and leads to a Maxwell term 
for gauge field a^,, due to the mutual CS coupling. Thus 
the spinon Za is coupled to a dynamical gapless U(l) 
gauge field Ufi. Then the GSM of the condensate of Za is 
S^/U{1) — S"^, which is equivalent to the collinear Neel 
order. Thus under the assumption s„ < 0, Eq. [5] reduces 
to the CP(1) model that describes the deconfincd QCP 
between Neel and VBS orde r^'^i^^ . 

We have shown that the mutual CS formalism Eq. [S] 
unifies many previously discussed exotic states and ex- 
otic phase transitions. A more detailed discussion of the 
phase diagram can be found in Refi^. 



B. Z2— Neel transition driven by (e,m) excitation 

Now we are ready to discuss our theory for the direct 
continuous transition between Z2 liquid phase and Neel 
phase. In a Z2 topological phase, using the standard 
notation, there are three types of topological excitations: 
the electric excitation e, the magnetic excitation m, and 
their bound state (e, m). In Eq. |6l the spinon field Za 
is the e— type excitation, while the vison field v is the 
m— type excitation. Eq. [S] is based on the assumption 
that inside the Z2 liquid phase the e— type and to— type 
excitations have lower energy than {e,m), thus in the 
global phase diagram Fig. [TJ the Neel and Z2 topological 
phases are separated by a multicritical point Sz = = 0. 
However, if we consider the opposite possibility, namely 
the (e, m)— type excitation has the lowest energy in the 
Z2 spin liquid, then a different quantum phase transition 
can occur by condensing the (e, m)— type excitation. 

Let us first take the simplest Toric code modeli^ 
as an example: H = E» -<-:.<:.<-y<y - 
'^i,x'^i,v'^i+x,y'^i+y,x- The "condensation" of an excita- 
tion simply means that the system enters a phase were 
the kinetic energy of this excitation dominates. It is well- 
known that in the Toric code model the condensation 
of the e— excitation is driven by a magnetic field hzaf ^, 
while the condensation of to— excitation is driven by field 
hxcrfi^, because these two fields enable the hopping of e 
and TO excitations respectively. In order to "condense" 
the (e, m) excitation, we simply need to turn on field 
hyof ^, which hops the (e,TO) excitation along the diago- 
nal directions of the square lattice. When any of the three 
excitations is condensed, the system enters a trivial po- 
larized state without any topological degeneracy. Gener- 
ally speaking, in the topological phase, starting from one 
of the topological sectors on the torus, the other sectors 
can be generated by locally creating a pair of topologi- 
cal excitations, and annihilating them after adiabatically 
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FIG. 2: (a) For a SU(2) invariant system, if one condenses the 
ordinary spinon of the Z2 spin liquid phase, the spin ordered 
state must have ground state manifold SO(3). One example 
state of this kind is the spiral spin density wave. (6) If the 
(e, m) type of excitation of the Z2 spin liquid condenses, the 
spin order can be the ordinary Neel order. At the ^2— Neel 
QCP, both Neel and columnar VBS order parameters have 
power-law correlation. 

moving one excitation of the pair around the torus. Be- 
cause all three types of topological excitations are mutual 
semions, condensing one of the three excitations will lead 
to a strong local flux fluctuation for the other two exci- 
tations, thus the other two excitations are confined in 
this condensate, i.e. the system no longer has topological 
degeneracy. 

In our current case, both e and m excitations carry 
extra global symmetries besides their gauge charges. In 
order to describe the (e, to) excitation in our situation, 
let us define new complex bosonic fields cjja and ipa'- 

(f>a=ZaV, i>a=ZaV*. (8) 

(j)a and ipa carry the quantum number of (e, to) excita- 
tion. Because w is a complex variable, fields 0^ and 4'a 
are independent from each other, and they interact with 
each other as follows: 

C = \{df, - iUf^ - ibf,)(j>a\'^ + \{df^- iaf_, + z6^)-0a|^ 

a 

+ r{\(pa\'^ + \ipa\'^) + -ef,,,paf,d^bp 

TT 

Notice that (jja and ipa carry gauge charges of both gauge 
fields and b^. In order to understand the QCP at r = 
more quantitatively, it is more convenient to define new 
gauge field = a,, ± 6^, then the Lagrangian reads : 

^ = J2\{d,-iA;M' + \{d,-zA-)ij^\' 

a 
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+ 9{W f +g{H\'f +u\(I)\^\tP\'' ~w^'^a(l)-Tp^di>. (10) 

In this field theory, (pa and V'a are almost decoupled 
from each other, i.e. they are only coupled through the 
quartic terms u and w. The mass gaps r for (pa and tpa 
are equal, because the vison modes v and v* are guar- 
anteed to be degenerate by the symmetry of the square 
lattice. (j)a and ipa are introduced as bosonic fields, but 
gauge fields v4+ and A~ make them fermionic fields af- 
ter the standard flux attachment, due to the existence of 
the Chern-Simons terms in this Lagrangian. In our for- 
mulation, fields (/)q and ^a can still condense by tuning 
parameter r in Eq. 1101 After (pa and ipa both condense 
simultaneously, ^+ and A~ are both Higgsed, and in the 
Higgs phase the only gauge invariant operators are 

(l)^ap, ip'^ai). (11) 

Since these two vectors both carry the same quantum 
number as the Neel order parameter z'^az, in Eq. [TU] w 
is naturally positive, thus these two vectors are aligned 
parallel with each other, so the condensate of (pa and tpa 
has a manifold S"^ , i.e. it is the standard collinear Neel 
order. The difference between this new transition and 
the ordinary spinon theory is illustrated in Fig. [2] 

What is the universality class of this transition? The 
simplest possibility is that, both u and w are irrelevant at 
the transition, although they are relevant in the conden- 
sate of pa and tpa- Under this assumption (pa and ipa are 
completely decoupled at the transition r = 0. Then in 
this case this transition is described by the simple Chern- 
Simons-Higgs model: 

N 

iN 

+ -^^iivpA^_idyAp. (12) 

Here we have generalized the equation to have N flavors 
of matter flelds pai and introduced a statistical angle 
0. Our physical situation corresponds to iV = 2 and 
9 = TT. Notice that Eq. [12] explicitly breaks the time- 
reversal symmetry due to the Chern-Simons term. But 
the complete theory Eq.fTUlis time-reversal invariant, be- 
cause under time-reversal transformation (pa and ipa are 
exchanged, the two gauge flelds v4+ and A~ are also ex- 
changed. 

The critical exponents of this transition can be com- 
puted using a systematic 1 /N expansion. Ref has com- 
puted the critical exponent v defined as ^ ~ here 
we will focus on the scaling dimension of (p'^T°'p at the 
QCP, where T"" is the SU(N) generator. To the first order 
1 /N expansion, this scaling dimension reads 

A[^tT^^] = l + Af---^^V (13) 



Let us briefly comment how we obtain this result. Similar 
calculation without the 9 term was obtained before. See 
Fig. 3 and Fig. 4 of the previous work^ - for necessary Feyn- 
man diagrams. First, we need to evaluate wave function 
renormalization of p from both gauge fluctuation and 
the density fluctuation, which contain the factor(l/A^). 
Then, using the standard operator insertion method, 
one can calculate renoramlization function of the corre- 
sponding vertex. We note that the traceless condition, 
(Tr(T") = 0), reduces one diagram compared with the 
calculation of the scaling dimension of and simplify 
our calculation. 

In the limit of 6* ^ oo, i.e. the CS term is ef- 
fectively zero, these results converges to the ordinary 
CP(iV — 1) results computed in Ref.?^. In the limit of 
6* 0, the gauge fluctuation is totally frozen by the CS 
term, and the universality class of this quantum criti- 
cal point only acquires corrections from the short range 
self-interaction between field pa^ thus it is equivalent to 
an 0(2Af) transition of the 0(2A^) bosonic vector field 
(Re[0i], • • • Re[(/)Ar], lm[0i], ■ • • Im[(/)Ar]). Scaling dimen- 
sion of the "Neel" type operator (A[0^^r°0]) in our theory 
is larger than that in the CP(A^ — 1) theory with large 
TV, i.e. at the Z2— Neel QCP, the anomalous dimension of 
the Neel order parameter is predicted to be larger than 
that of the deconfined QCP between the Neel and VBS 
order. This prediction can be tested in the future by a 
careful comparison between the critical exponents of the 
Ji — J 2 model and the J — Q modeP^"-^. 

It is pretty clear that at least in the large— limit, 
the perturbation of u in Eq. [TO] is irrelevant, because in 
this limit the scaling dimension A[|(/)p] = A[|7/'p] = 2, 
i.e. A[u] = —1. Higher order 1/7V or e expansion is de- 
manded to determine whether w is relevant or not at this 
transition. 

Assuming at the QCP r = both u and w are irrele- 
vant, then besides the Neel order parameter, some other 
physical order parameters also have power-law correla- 
tion. For example, the columnar VBS order parameter 
can be written as 

VBS ^ ^IPaMa «'A^a, (W) 

where Ma is the monopole operator for gauge field a^. 
When pa and 'ipa both have a large A'' component, the 
scaling dimension of A^a is proportional with A. Thus 
with large A^ the VBS order parameter is expected to 
have a much larger scaling dimension compared with the 
Neel order parameter at the Z2— Neel QCP. We stress 
that the VBS order parameter has short-range correla- 
tion in the Z2 spin liquid and the Neel phase, its emergent 
quasi long range correlation occurs only at the QCP. This 
result has already been confirmed numerically in Refi^, 
and it was demonstrated that the scaling dimension of 
the VBS order parameter is indeed larger than that of 
Neel order at the QCP"*. 

In 2-1-1 dimension, the entanglement entropy of a con- 
formal field theory can in general be written as S' = 
cL — pi, where the first term is the nonuniversal area law 
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contribution, while the second term is a universal con- 
stant. In Ref<^, it was argued that at a QCP where 
a bosonic field condenses while coupling to a discrete 
gauge field, the universal entanglement entropy is a di- 
rect sum of the contribution from the bosonic matter 
field and the contribution from the discrete gauge field: 
/3 = /3f, -f /3 gauge- This conclusion is based on the as- 
sumption that the matter field dynamics is not affected 
by the discrete gauge field in the infrared limit, and this 
is indeed true for the XY* transition observed in Ref^iS. 
However, at the exotic Z2— Neel transition discussed here 
where the (e, m)— type excitations condense, the bosonic 
matter fields 0^ and ■0a are indeed strongly affected by 
the gauge field, thus at this QCP the universal entan- 
glement entropy /3 is no longer a direct sum of the two 
different degrees of freedom of the system. The univer- 
sal entanglement entropy of field theory Eq. [12] in the 
large— limit can be found in Refi^. 

C. A Toy model with iV = 1 

Now let us discuss a toy model with iV = 1. This is 
actually the case where the critical theories can be all 
understood exactly. This field theory with iV = 1 can be 
applied to the following extended Toric-code model: 

i 

Here the e— type (m— type) excitation is the end of a 
string product of cr^ (c^). The e and m— type excita- 
tions view cr^ and as Z2 gauge fields respectively, and 
the hx and terms enable the hopping of these excita- 
tions. Unlike the standard toric-code modeU^, here we 
keep Kx, > 0. When K^, > 0, both cr^ and a"^ 
have a tt— flux in the ground state. Then the dynamics 
of both e and m type of excitations are frustrated, and 
both excitations have two different minima ±Q in their 
band structure. As a result, the low energy dynamics of e 
and m excitations are described by complex scalar fields 
z and V expanded at momentum Q. If one of these two 
fields condenses while the other one remains gapped, the 
Z2 topological order is destroyed, and the system must 
spontaneously break the lattice translation symmetry as 
well. The condensates of e and m— type excitations phys- 
ically corresponds to the valence bond solid phase of cr^ 
and (7^ respectively. 

When = 1, if e and m type of excitations condense 
separately, then the phases in Fig. [T] would be Z2 liquid, 
VBS order of a^, trivial phase, and VBS order of cr^ 
(counted counterclockwise around the multicritical point 
Sz = Sv = 0). On the other hand, if the bound state 
(e, to) has the lowest energy in the Z2 liquid phase, then 
again we can introduce two independent complex fields 
(j) and Tp as 4> = zv, ip = zv* . Then the transition driven 



by the condensation of (j> and ■0 is described by Eq. [12] 
with iV 1 and 9 tt/2. 

What kind of transition is this? If in Eq. [T^] the com- 
plex field is also coupled to an external U(l) gauge 
field A'jf*, then we can see that in the disordered phase 
of 0, after integrating out the massive and dynam- 
ical gauge field Af^, the lowest order contribution to 
the effective Lagrangian of A^f^ is still a Maxwell term: 
Ceff ~ ((9A^^*)2 + c6'{92A«^*)(a^e^t) + . . .. While in 
the condensate of 0, the effective Lagrangian of A'^f^ ac- 
quires a Chern-Simons term at level 1. This analysis 
implies that this transition is equivalent to a topological 
transition between a trivial insulator and a Chern insu- 
lator with Chern number 1 . The universality class of this 
type of topological transition of Chern insulator is very 
well-understood, it can be simply described by a 2-t-ld 
Dirac fermion: 

C = i)^^df_,4> -f TO'0'0, (16) 

here the trivial insulator and Chern insulator correspond 
to TO > and m < respectively, and m — Q corresponds 
to the quantum critical point at r = in Eq.[T2]with N = 
1. Thus we conjecture that when TV = 1 and = it/2, the 
critical point in Eq. [12] is dual to a massless free Dirac 
fermion. In Ref a similar conjecture was made that 
the 3D XY transition is dual to a massless Dirac fermion 
coupled to a noncompact U(l) gauge field. 



IV. SUMMARY AND DISCUSSION 

In this work we have discussed a possible theory for the 
direct continuous transition between the Z2 liquid phase 
and the Neel order, and this is a candidate theory for the 
liquid-Neel transition observed in Rcf.'^'^^ . We have taken 
the square lattice as an example, but results discussed in 
this paper can also be applied to the honeycomb lattice 
after straightforward generalization. 

In our theory, we exploited the fact that in two di- 
mension, the e and m— type excitations are both point 
like defects, thus their nontrivial statistics can be de- 
scribed well with a mutual Chern-Simons theory. By 
contrast, in a three dimensional Z2 liquid phase, there is 
a mutual semion statistics between the point particle like 
e— excitation and loop like to— type excitation. Thus the 
effective field theory for the three dimensional Z2 liquid 
phase is the so-called BE theory Ceff ^ ■^^fj.upTafj.d^bpr, 
where is the U(l) gauge field that couples to the 
e— type point particle, and is an antisymmetric rank- 
2 antisymmetric tensor gauge field that couples to the 
TO— type loop excitation. The global phase diagram 
around the three dimensional Z2 liquid phase is another 
interesting subject, and we will leave it to future studies. 
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